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1. Introduction 



In recent years there has been a resurgence of interest in charmonium physics, both theoreti- 
cally and experimentally. The charmonium states below threshold are considered important topics 
for study by lattice QCD. They provide us with a testing ground for lattice heavy quark methods 
which are also used to determine elements of the CKM matrix involving b quarks. New experi- 
ments at CLEO-c aim to confront lattice QCD calculations with extremely precise experimental 
data from the charm sector. In addition, the discovery of numerous new, and unexplained, bound 
states such as the X(3871) and the presence of theoretically allowed exotic states has renewed 
efforts in charm physics. 

Lattice QCD can in principle answer many of the outstanding questions but to do so requires 
high precision numerical simulations. 

In these proceedings we present preliminary results from a relativistic simulation of the char- 
monium spectrum, using all-to-all propagators on a dynamical anisotropic lattice. An earlier study 
of the spectrum on a smaller volume was presented in Ref. [jl|]. The charmonium spectrum at finite 
temperature was also determined in this way and is presented in Refs [^, ||. |J] . 



2. The dynamical anisotropic action 

The gauge action is a two-plaquette Symanzik-improved action. It is written 



S - P 



where Q.ff refer to simple and rectangular plaquettes in the spatial and temporal directions and %® 
is the bare anisotropy. The full details are described in Ref [g]. 
The anisotropic fermion action [||] is written 
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where a s and a, (a s ^> a,) are the spatial and temporal lattice spacings respectively and ^ is the 
bare anisotropy. 

In an anisotropic simulation the bare anisotropics, which are inputs in the gauge and fermion 
actions must be tuned so that the measured value of the anisotropy in a simulation takes its "target" 
value. In the quenched theory the gauge and fermion anisotropics can be separately tuned: however 
in the dynamical theory these must be simultaneously tuned. The details of this tuning for quenched 
and dynamical QCD are described in Refs. [||, 



3. Simulation details 

We performed simulations at Nf = 2 on N s x N t = 12 3 x 80 lattices with 250 configurations at 
a sea quark mass close to the strange quark mass. The simulation details are summarized in Table 
The target (renormalised) anisotropy in this study is 6. We found that, in contrast to the quenched 
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Configurations 250 (a t m c = 0.117,c? f m sea = ^mug^t = —0.057) 



Dilution 

Physics 

Volume 

N f 
a s 

mx/nip 



time and space even/odd 

S,P,D waves, S wave radial excitations, and hybrids 

12 3 x 80 

2 

~0.17fm 

-7GeV 

~0.55 



Table 1: Simulation details for this study of charmonium. 



case [||] the anisotropy in the charm sector had to be tuned separately to that in the light quark 
sector to obtain the same renormalised anisotropy. This mass-dependence of E, in the dynamical 
theory may be a discretisation effect and is under investigation. The charm quark mass used in the 
simulation, a t m c = 0.117, was tuned to obtain the correct JpV mass. We use all-to-all propagators 
with dilution, no eigenvectors and two noise vectors, as described in Ref. [Q], for better signal to 
noise ratios. An additional advantage of all-to-all propagators is that constructing an extended basis 
of operators for better overlap with the states of interest is easier than with point propagators. The 
list of operators used in this study is given in Table ^[ The spin-averaged (IP-IS) splitting is used 
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Table 2: The basis of operators. The notation used for the gluonic paths, denoted with s,-,i = 1,2,3, is 
defined in Ref. ©. 



to set the lattice spacing, yielding a t = 0.028fm and a s = 0. 17fm. 

Since time dilution introduces a random noise source at each timeslice, effective masses fluctu- 
ate more from timeslice to timeslice than is usually seen with point propagators making it difficult 
to identify a good plateau region. However, these fluctuations do not affect exponential fits and 
can be reduced further by increasing the dilution level. A more accurate reflection of the quality 
of the data and stability of the fits can be seen from a sliding window plot. For a fixed value of 
?max the fitting window, (? m in,fmax) is varied by changing t m \v The sliding window plots show the 
fitted masses at each interval as a function of t m \n . Sliding window plots for IS and 2S states of 
J/\ff are shown in Figure |]. The plots show that the IS fitted mass is stable as fmin is varied over 
twenty timeslices and the 2S mass is stable over seven values of t m \ n . For all the fitted masses 
shown here the quality of the fit is determined and the stable regions also have good % 2 /N±f_. The 
other charmonium states studied also show extremely stable fitted masses. Based on this analysis 
the hyperfine splitting is approximately 30 MeV in contrast to the experimental value of 1 17MeV. 
The reason for this underestimate is not known yet and will be investigated in future work. In par- 
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Figure 1: Sliding window plots for the ground and first excited states of the J/\j/ at zero momentum. The 
solid lines show the final best-fitted masses chosen. 



ticular, the lattice spacing in this study is quite coarse and the action does not include a £ • B term. 
In addition it is known from previous work that the hyperfine splitting increases with decreasing 
lattice spacing and for lighter sea quark masses in dynamical simulations. These improvements 
will be included in future simulations. 

To determine the radially excited states we use a variational basis of operators included ex- 
tended operators. Our results for the ground state mass using this approach agrees within two 
percent with the same state determined from exponential fits. Two smearings and two different 
operators were used to increase the variational basis when determining the r] c and The first 
excited state of the r] c and J/W are m^s) = 3645 ±12 MeV and m^^s) = 3701 ± 12 MeV re- 
spectively. This analysis is being extended to the P and D waves. Our preliminary charmonium 
spectrum, including the P- waves, D-waves and the 1~ + hybrid is shown in Figure ||. 



4. Disconnected diagrams 

In most lattice calculations of the charmonium spectrum and the hyperfine splitting in par- 
ticular, the effects of disconnected diagrams in the two point correlator are ignored. It has been 
estimated that these OZI-suppressed diagrams may in fact contribute ~ 20 MeV to the hyperfine 



splitting jjllj, |12J], bringing traditionally low lattice determinations into better agreement with ex- 
periment. However, the inclusion of such diagrams is not easy, requiring all-to-all propagators 
which in the past have resulted in noisy data from which it is difficult to extract a convincing sig- 
nal. It is hoped that the all-to-all propagator algorithm used here may solve this problem using 
the appropriate dilution to reduce noise. In this first attempt the diagrams were directly calculated 
using the same dilution already described (time and space). We found that the disconnected dia- 
grams are still noisy and in fits to the full correlator, given by the difference of the connected and 
disconnected correlators the fit window is reduced to ten data points. Sliding window plots for 
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Figure 2: The cc spectrum which is normalised to the n c mass. The scale has been set from the (IP — 15) 
splitting. Experimental values, where known, are denoted with right-facing triangles. 



the full correlators, Cf u n(7) = Connected (0 — Qisconnected(O) we shown in figure |3[ The plots show 
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Figure 3: Sliding window plots showing the fitted mass in lattice units for the full correlators, Ct u u = 
Cconn — Cdis for the tj £ and J '/ip. The plots show the effect of changing f max in the fits. The poor quality 
of the signal means f ma x has been moved in from the centre of the lattice approximately timeslice ten. The 
straight line represents the ground state mass obtained from the connected parts of the correlators. 



no discernible difference after including the disconnected contributions, however we are currently 
investigating other dilution schemes to optimise the signal. 

5. Conclusions 

We have presented our preliminary results from a simulation of the charmonium spectrum on 
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Nf = 2 dynamical 12 3 x 80 anisotropic lattices. All-to-all propagators and a variational basis of 
operators have been used the determine the spectrum of S, P and D waves, their radial excitations as 
well as the 1~ + hybrid with good statistical precision. In addition, we have investigated the effects 
of the disconnected diagrams on the hyperfine splitting. We find a negligible difference. However, 
these preliminary results need to be investigated further. In particular we are planning to explore 
dilution strategies and study the volume dependence of the higher excited states. 
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